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Abstract 

The advent of increasingly precise gyroscopes has played a key role in the technological develop- 
ment of navigation systems. Ring-laser and fibre-optic gyroscopes, for example, are widely used in 
modern inertial guidance systems and rely on the interference of unentangled photons to measure 
mechanical rotation. The sensitivity of these devices scales with the number of particles used as 
1/s/N. Here we demonstrate how, by using sources of entangled particles, it is possible to do bet- 
ter and even achieve the ultimate limit allowed by quantum mechanics where the precision scales 
as 1/N. We propose a gyroscope scheme that uses ultra-cold atoms trapped in an optical ring 
potential. 

PACS numbers: 03.75.-b, 03.75.Dg, 03.75.Lm, 37.25. +k, 67.85.-d 



I. INTRODUCTION 



Optical interferometers have revolutionised the field of metrology, enabling path length 
differences to be measured, for the first time, to less than the wavelength of the light being 
used. As well as their more familiar linear versions, interferometers can also be used in ring 
geometries to make accurate measurements of angular momentum. These interferometric 
gyroscopes surpass the precision of their mechanical counterparts and form a key component 
of many modern navigation systems. They work by exploiting the different path lengths 
experienced by light as it propagates in opposite directions around a rotating ring. For 
instance, in the Sagnac geometry p[| photons are put into a superposition of travelling in 
opposite directions around a ring and, when rotated, the two directions acquire different 
phases. This phase difference is directly related to the rate of rotation and can be measured 
by recombining the two components at a beam splitter and recording the intensity at each 
of the outputs. 

Such schemes use streams of photons that are independent of one another, i.e. not 
entangled. In this case the measurement accuracy is fundamentally limited by the jitter 
in the recorded intensities due to the fact that photons come in discrete packages. This 
is known as the shot-noise limit and restricts the measurement to a precision that scales 
inversely with the square root of the total number of photons. One possible way to beat 
this precision limit is to use entangled particles [2J. In fact it has long been known that the 
precision of optical interferometers is improved with the use of squeezed or entangled states 
of light pHH]. In principle, this should enable us to reach the Heisenberg limit whereby the 
precision scales inversely with the total number of particles. 

The precision of gyroscopic devices can be improved further still with the use of entan- 
gled atomic states due to their mass enhancement factor over equivalent photonic devices 
[S]. The use of entangled atomic states to make precision measurements of rotations was 
first proposed in [6] . Since then, and with the experimental realization of Bose-Einstein con- 
densates, the use of entangled atoms for precision measurements has been widely researched. 
As in the optical case, of key importance to the ultimate precision afforded by an atomic 
device is the input state used. The use of number-squeezed atomic states will allow for 
Heisenberg limited precisions and as such there has been much research into the generation 
and uses of these squeezed states J7HTU]. Several proposals have already been made to use 
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these squeezed, and other entangled atomic states, such as maximally entangled 'NOON' 
states, to make general phase measurements with sub-shot noise sensitivities [TTHTT] . It 
has also been shown that uncorrelated atoms can also achieve sub shot-noise sensitivities of 
rotational phase shifts [18J by using a chain of matter wave interferometers, or a chain of 
gyroscopes. 

Here we propose a gyroscope scheme that uses squeezed and entangled atomic inputs to 
push the sensitivity of rotational phase measurements below the shot noise limit. We extend 
the investigation of optimal input states to an experimentally accessible atomic gyroscope 
scheme capable of measuring small rotations. It works by trapping ultra-cold atoms in a one- 
dimensional optical lattice in a ring geometry and carefully evolving the trapping potential. 
We investigate the precision achieved with different inputs and show that although a so- 
called NOON state produces the best precision, another state (sometimes referred to as a 
'bat' state) that is created by passing a number-squeezed state through a beam splitter, is 
far more robust and might therefore be a preferred candidate. 

It should be noted that it is possible to beat the Heisenberg precision scaling of 1/N 
in some measurement schemes. In fact precision scalings of iV _3//2 can be achieved even 
when the initial state is unentangled in a few specific metrology protocols [191 ED]- This is 
achieved through a nonlinear coupling between the quantum probe and the parameter to be 
measured. Therefore, BECs with their particle interactions, may naturally lend themselves 
to this. However, the challenge would be to find ways of coupling the angular momentum we 
wish to measure to the scattering length of the atoms. While this may provide an interesting 
future direction to this work, in light of this difficulty, we concentrate for now on attaining 
the Heisenberg limit through optimizing the input state. 

II. THE SYSTEM 

Our system consists of a collection of ultra-cold atoms trapped by the dipole force in an 
optical lattice loop of three sites. Rings with this geometry have already been experimentally 
demonstrated [2H [22]- For a sufficiently cold system, we need only consider a single level in 
each site and so can describe the system using the Bose-Hubbard Hamiltonian, 



H 



2 



2 2 



Yl J i • • + a j- a j) + Yl y 3 a f a ) 



(i) 



j=0 




3 



where cij is the annihilation operator for an atom at site j and Jj is the coupling strength 
between sites j and j + 1. The ring geometry means that aj = 0^+3. The parameter Vj is 
the strength of the interaction between atoms on site j and ej accounts for the energy offset 
of site j. In general, we take the zero point energy to be the same for each site and so set 
€j = 0. For the purposes of this work it will be convenient to describe the system in terms 
of the quasi-momentum (or flow) basis which is related to the site basis by, 

«, = ^X> Wfl - ( 2 ) 

where a%. corresponds to the annihilation of an atom with k quanta of flow. We shall use 
positive subscripts to refer to clockwise flow and negative subscripts to refer to anti-clockwise 
flow. 



III. SCHEME 1: UNCORRELATED PARTICLES 



The first scheme we present uses unentangled atoms to achieve shot noise limited pre- 
cisions. To begin with, the potential barriers between the sites are high and N atoms 
are contained within one site, say site zero. The initial state of the system is therefore 
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\N, 0, 0) where the terms in the ket represent the number of atoms in sites zero, one 
and two respectively. 

The first step is to rapidly reduce the potential barrier between just two sites, we choose 
sites zero and one, in such a way that the two sites remain separate but there is strong 
coupling between them. This must be done rapidly with respect to the tunneling time, 
but slowly with respect to the energies associated with excited states in order to ensure 
the system remains in the ground state. This separation of timescales has already been 
demonstrated experimentally |23j. In this regime the coupling between the two sites is 
much larger than their on site interactions and the Hamiltonian describing the two sites is, 
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Importantly, the remaining two barriers are high (V ^> J) and so prevent tunneling between 
sites one and two and sites two and zero. 

The system is left to evolve for time t = n/AJ whilst this barrier is low. This is equivalent 
to applying a two port 50:50 beam splitter to our initial state (as shown in |15j) and so 



transforms \ip)uo to, 

Wia = -^=(aS + wOT0 > 0,0>- ( 4 ) 
Each individual atom is now equally likely to be on site zero or one. In other words we have 
N single-particle superpositions on the two sites. 

The next step is to apply a three port beam splitter, or tritter. This splitting procedure 
is described in detail in [21]. Essentially to achieve a tritter in this system we immediately 
lower the two remaining potential barriers, on the same timescale as before, and allow the 
system to evolve for a further t = 2tt/9J. This tritter operation is given by, 
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from which we can see \ip)ui is transformed to, 
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At this point we rapidly raise the potential barriers, 'freezing' the atoms in the lattice sites. 
Comparing \ip)u2 with equation §2§ we see that applying a 2tt/3 phase to site two results in 
a superposition of the a_i and a\ flow states. This phase is achieved by applying an energy 
offset, €2, to site two, whilst the barriers are high, for time t t = 47r/3e2- Offset application 
times of 500ns have been demonstrated experimentally [25] and it is this time we shall use 
in section VI C| when we assess the impact of non-zero interactions. 



We then immediately lower the barriers again so the atoms can flow around the loop. 
The resulting superposition can be written as, 

Win = yp^( a -i + °> °> ( ? ) 

where the terms in the ket now represent the number of atoms in each of the possible flow 
states, a and aii respectively. This is now a iV single particle flow superposition. 

At this point the a_i and a% states are degenerate, so |V ; )i/3 does not evolve. However, 
we now apply the rotation we wish to measure, u, to the ring which causes a phase, 9, 
to be applied around it. The energies of the two flow states now change according to the 
Hamiltonian, 
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After a time t w , and ignoring global phases, the state has evolved to 

and so a phase difference of <p = 2\/3Jt UJ sin(0/3) is established between the two flows. 

We now wish to read-out this phase difference from which we can directly determine u 
since u = h8/(L 2 m) where m is the mass of the atom and L is the circumference of the 
ring. The read-out procedure involves sequentially undoing all the operations performed 
prior to the phase shift. This is analogous to standard Mach-Zehnder interferometry where 
an (inverse) beam splitter is placed after the phase shift to undo the initial beam splitting 
operation. 

The undoing process begins with the application of a — 27r/3 phase to site two giving, 

1^5 = ((4 + e^ 3 a\ + e^al) + ie»{e*"<*4 + a\ + e^aj))^ |0, 0, 0) 

(10) 

where the terms in the kets now, once again, represent the number of atoms in sites zero, 
one and two. Next we undo the fritter by applying an inverse fritter, S^ 1 = Sg. The inverse 
fritter operation is discussed in detail in j21] but essentially is achieved by lowering all three 
barriers and allowing the system to evolve for time t = Air/QJ (i.e. twice as long as for a 
fritter) giving, 

W)m = (4 + ^4)^10,0,0) (11) 

which is equivalent to \ip)ui but with a phase difference, <fi. 

Finally we apply an inverse two port 50:50 beam splitter. This is achieved in just the 
same way as the two port beam splitting operation described above, but with a hold time 
of t = 37r/4J rather than t = n/AJ. The resulting state is, 

= -j= (cos (?\ (4) - sin (£j (a\)j |0, 0, 0) (12) 
meaning the probabilities of detecting each atom at site zero and site one are, 

Po = cob 2 (|) (13) 

Since the atoms are independent the total number detected in the two sites is given by a 
binomial distribution. The mean number of atoms detected at site zero is therefore (n ) = 
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iVcos 2 (0/2) and at site one it is (ni) = A^sin 2 (0/2). By counting the number of atoms 
detected at each site we can determine 0, and hence u, just as in a typical Mach-Zehnder 
interferometer. 

The precision with which this scheme enables us to measure u can be found by calculating 
the quantum Fisher information, Fq. This is a tool for evaluating the precision limits of 
quantum measurements and is independent of the measurement procedure. For a pure state 
|^(0)) it is given by [26] . 

_(* / (0)|* / (0)>-|(* / (0)|*(0))| 2 J. (14) 

We convert this into an uncertainty in <ft using the Cramer- Rao lower bound [27H29] , 

A<f>>l/y%. (15) 

Using \ip)u4: we find the maximum resolution scaling of <\> with N is N~ 1 ^ 2 , or equivalently 
A6 ~ y/3/(2Jt LJ cos(6/3)y/N). This translates to an uncertainty in uj of, 

Au; ~ ( (16) 
\L 2 mJ 2Jt w \fN 



where we have made the approximation that 6/3 <C 1. This has the well-known 1/yN 
scaling that is a signature of the shot-noise limit. 

To summarize, shot-noise limited measurements of rotations are made as follows: 

1. Apply a two port 50:50 beam splitter to the first two modes of the state \N, 0, 0). 

2. Perform a three port beam splitter (fritter) operation to the state. 

3. Apply a 27r/3 phase to site two. 

4. Leave the system to evolve for time t u under the rotation, u. 

5. Apply a — 27r/3 phase to site two. 

6. Perform an inverse fritter operation on the state. 

7. Apply an inverse two port 50:50 beam splitter to the first two modes. 

8. Count the number of atoms in each site. 

We will now show how our scheme can be modified to create entangled states and will 
investigate the effect of this entanglement on the precision scaling of our rotation measure- 
ments. 
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IV. SCHEME 2: THE BAT STATE 



We begin with N/2 atoms on site zero and on site one, i.e. |V>)bo — |iV/2, iV/2, 0). The 
production of dual Fock state BECs in a double well potential was first proposed in [30J. 
This number squeezed state could be achieved by slowly applying a double well trapping 
potential to a condensate so that a phase transition occurs to the Mott insulator state and 
has been demonstrated experimentally in three-dimensional optical lattices [5T] . Initially 
the barriers between the three sites are high and, as in scheme 1, a two port beam splitter 
is applied to the first two modes of \iP)bo- The resulting output is sometimes referred to as 
a 'bat' state since a plot of the amplitudes in the number basis resemble the ears of a bat. 

Steps two to seven are identical to scheme 1 resulting in, 

^ BA = VW(N/2)l ( (a " l)2 + Vl)2 ) |0 ' °' 0) (17) 

and 

1/2 2\ N / 2 

W™ = 2*(N/2)\ (( a °~ iQ + e** ("»"4 + aJ) J |0,0,0> (18) 

where a global phase has been ignored and <fi is again given by (j) = 2\/?>Jt U j sin(d9/3). This is 
very similar to the scheme in [16] , where a bat state is used to measure a phase difference in 
a general Mach-Zehnder interferometer set-up, and in fact results in the same output. The 
difference between the schemes is that ours has been adapted to measure a rotation around 
a ring of lattice sites rather than a general phase between two paths. 

To determine (f) we could count the number of atoms detected at each site and repeat. 
However, this requires nearly perfect detector efficiencies [32J and so is experimentally chal- 
lenging. Instead we use the read-out scheme described (in detail) in [16] to make our 
measurements. 

Essentially, after step seven the system is then left to evolve with the barriers high for r = 
7r/16V (note in the original paper r = 7r/8U because U = 2V here). The trapping potentials 
are then switched off and after some expansion time interference fringes are recorded. These 
two steps are the read-out steps and are what we shall refer to as step eight. The scheme is 
repeated many times and the visibility, V, of the fringes is calculated as in [16]. From these 
visibility measurements we determine <ft directly. Note that the quantum Fisher information 
depends only on the final state of the system and not on the measurement procedure. 
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As such we use \iP)ba to determine the precision scaling of 9 with N and find A6 



v/3/(2Jt w cos(9/3)^N(N/2 + 1)). This means the uncertainty in cu is, 

V^m/ 2Jt w cos(0/3)y/N(N/2 + l) \L 2 m ) V^Jt^N' 
where we have made the approximations N ^> 1 and 6/3 <C 1. This has the same number 
scaling as the Heisenberg limit. 

We now present a third scheme which offers a slight improvement in the precision scaling 
of our measurements. We then consider the advantages and disadvantages of schemes 2 and 
3 and discuss their experimental limitations. 



V. SCHEME 3: THE NOON STATE 



This scheme is very similar to scheme 1, the only difference is that the two port 50:50 
beam splitter (and its inverse) is replaced with a two port quantum beam splitter (and 
an inverse two port quantum beam splitter). A two port quantum beam splitter (QBS) is 
defined device [33] that outputs (|JV,0) + e^|0, N))/y/2 when |JV, 0) is inputted. The 
terms in the kets represent the number of particles in the two modes. This output is called 
a NOON state. 

As in scheme 1, the three potential barriers are initially high and N atoms are inputted 
into site zero, \iP)no — I A, 0,0). The first step is to apply the QBS. This is done using a 
scheme proposed in [T5] . 

The QBS begins with the application of a two port 50:50 beam splitter between sites 
zero and one, as previously described. A 7r/2 phase is then applied to one of the two sites 
using an energy offset as above. At this stage K> J and the interactions are tuned such 
that their strength on one site is an integer multiple of the strength on the second site [46J . 
The system is left to evolve for t = it /2V in this regime after which a second two port beam 
splitter is applied. These steps output, 

I^ati -> ^= ((al) N + e*(a\) N ) |0, 0, 0). (20) 

where £ is some relative phase established by the splitting procedure. 

Here we have a superposition of all N atoms on site zero and all on site one. At the 
equivalent stage in scheme 1 we had iV single particle superpositions (see equation [5]) . It is 
this difference that is responsible for the improved precision. 
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Steps two to six are identical to those in scheme 1 giving, 



and 



1^4 = ^= ((aU) N + e% iN t>(al) N ) |0,0,0> 
mm = -tL= ((ai) N + e% iN *(a\) N ) |0, 0, 0). 



(22) 



(21) 



Finally we apply an inverse two port quantum beam splitter to the system. This is 
achieved by sequentially undoing all the steps of the QBS. So, 

1. Apply an inverse two port beam splitter. 

2. Raise the barriers and tune the interactions as before. 

3. Leave the system to evolve for t = ir/2V. 

4. Apply a — n/2 phase to the same site as before. 

5. Apply a second inverse two port beam splitter. 

6. Raise the barriers. 

This results in all N atoms detected at site zero or all at site one with respective proba- 
bilities, 



By repeating the scheme many times, each time recording the site on which all N 
atoms are detected, 0, and hence u, can be determined. Using the quantum Fisher 
information and the Cramer-Rao lower bound the maximum resolution of scheme 3 is 



We see the NOON state offers a slight improvement in resolution scaling over scheme 2. 
It has the same number scaling, but the numerical factor is y/2 better. However, this slight 
improvement in resolution may come at great experimental expense. We now investigate 
the effect of experimental limitations on schemes 2 and 3. 

VI. COMPARISONS AND PRACTICAL LIMITATIONS 




(23) 



A9 ~ y/3/(2 Jt u cos(#/3)A0 meaning, for 9/3 < 1, 




(24) 



Both schemes 2 and 3 allow Heisenberg limited precision measurements of small rotations 
with the precision capabilities of scheme 3 being marginally favorable. Our descriptions 
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thus far, however, have only considered the idealized case and we have neglected important 
physical processes that may limit the experimental feasibility of the schemes. We now 
reevaluate the schemes when these limitations are accounted for. 

A. Particle loss 

We need to account for the fact that atoms may be lost during the schemes. It is well 
known that NOON states undergoing particle loss decohere quickly and so soon lose their 
Heisenberg limited sensitivity. However, it has previously been shown that certain entangled 
Fock states, that allow sub-shot noise phase measurements, are much more robust to these 
losses |34j . Here we wish to see how resistant the bat state is to particle loss in comparison 
to the NOON state by determining the precisions afforded by both schemes in the presence 
of loss. 

We investigate the effects of loss using a procedure described in [351 EE] which models 
loss from the middle section (between the two beam splitters) of an ordinary two path 
interferometer using fictitious beam splitters whose transmissivities, r], are directly related 
to the rate of loss. The equivalent loss in our schemes is loss from the momentum modes 
during t w . It was shown in [331 EE] that the time at which loss occurs during this interval is 
irrelevant. Since losses are equally likely from both modes we consider equal loss rates, rj, 
from each momentum mode. 

To compare the effects of particle loss on our two schemes we determine Fq and A0 for 
each scheme for different loss rates, or different rj. As shown in [351 EE] Fq varies with rj as, 



where I is the total number of atoms lost, l ai is the number lost from the ati mode, Pi a i-i 
is the probability of each loss event and 

N-(l-l ai ) 



N I 



fq-^fq y] pica ,1-icc, i6, 



(0)1 



(25) 



1=0 l ai =0 



(0)) 
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(26) 



Here 




(27) 
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and (3 m is 



for the bat state and 



Pr, 



•7V/2 



2 N / 2 {m/2)\{N/2 -m/2) 

/3m — 
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1/^2 for m — 0,N 
for m ^ 0, N 
for the NOON state. Fq is found numerically using 

F Q = Tr[p{4>)A 2 } 

where A is the symmetric logarithmic derivative which is defined by, 

d P {4>) l, 



(29) 



It is given by, 



(A) 
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(30) 



(31) 



Aj + Xj L _ T Jy 

in the eigenbasis of p(0), where Ay are the eigenvalues of p((f>). Combining these equations 
gives 



Aj + A,- 



d P (4>) 



(32) 



where are the eigenvectors of p(0). 

Figure [T] shows how A0 varies with rj for the two schemes when iV = 10. As expected the 
NOON state achieves the best precision when i] = 1, or when there are no losses. However, 
as f] decreases the bat state soon becomes the favored scheme. The lower bound of the 
shaded area is the Heisenberg limit and the upper bound is the precision achievable when an 
uncorrelated, or classical, input state is used, as in scheme 1. We see the uncorrelated state 
soon outperforms the NOON state. However, the bat state outperforms the uncorrelated 
input for approximately half the loss rates shown. Since it is unlikely half the atoms would 
be lost in an experiment, the bat state, unlike the NOON state, appears to offer an experi- 
mentally feasible increase in precision over classical precision measurement experiments. In 
the remainder of the paper we therefore assess the impact of experimental limitations on 
just scheme 2. 
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FIG. 1: The uncertainty of <f> for different rates of loss, rj, for N = 10. The blue dashed line shows 
A(p for scheme 2 and the green solid line shows A(f> for scheme 3. The upper bound of the shaded 
region is the precision afforded by scheme 1 (the classical precision limit - black dashed-dotted 
line) and the lower bound shows the Heisenberg limit. Scheme 3 soon becomes less favorable than 
scheme 1, whilst scheme 2 is much more robust to losses. 

B. Variations in N between experimental runs 

Our scheme requires many repetitions of the gyroscope procedure in order to build up 
interference fringes from which <fi can be determined. We have assumed thus far that each 
run involves exactly iV atoms. However, in an experiment N is likely to fluctuate between 
runs. The effect of fluctuations of order y/N on the bat state input are discussed in \Uj\ . 
In that case, an ordinary two-path linear interferometer is used but the same results apply 
here. It is shown that while the interference fringe signal is degraded by these fluctuations, 
the approximate Heisenberg limited sensitivity of the scheme is not destroyed. As expected, 
the larger N, the smaller the fluctuation effects, which is good since we would ideally work 
in the limit of large N since this gives the best improvement in precision. 

C. Interactions 

So far we have considered only the idealized system setting V = (apart from in the 
detection step, step eight, where we require large interactions to minimize small coupling 
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effects), and J = in the low coupling regime. While interactions can be tuned to extremely 
small values using Feshbach resonances it is an unrealistic assumption to discount them 
altogether. Likewise it is unrealistic to completely neglect coupling effects in the low coupling 
regime. Here we consider the effect of non-zero interactions and non-zero coupling strengths 
in the low coupling regime on scheme 2. To determine experimental orders of magnitude for 
V and J we use the approximations [38] . 

Tr 2a 9 V(?El 

V ps ° R 33 

-i-(-(?)/I)(/I+(/I) 3 ) 

where a s is the scattering length, Vq the barrier height, Er the recoil energy, A the wavelength 
of the lattice light and D the transverse width of the lattice sites. Feshbach resonances 
can tune a s to values smaller than the Bohr radius for some BECs [39] and in the high 
coupling regime barrier heights of order Vq = 2Er have been demonstrated |40j. Using 
light of wavelength A = D = 10/zm and 87 Rb atoms, interactions can therefore be tuned 
to V ~ 10" 3 Hz and J ~ 10Hz in this regime. While in the low coupling regime, where 
Vq = 35Er [23J, we find V ~ lCT 2 Hz and J ~ 10~ 2 Hz. Note that in the detection process 
the system evolves for t = ix/lQV with high potential barriers. Here we require large 
interactions to minimize small coupling effects. Taking a s = 9000ao [2] gives V ~ 100Hz. 
Using these values we assess the impact of non-zero interactions and coupling strengths. 

As iV increases the occupation number per site increases and as such the effect of non-zero 
interactions become more pronounced. We would therefore like to determine the maximum 
number of atoms our system can tolerate before these effects become too destructive. To 
do this we measure the fidelity between the output of the gyroscope in the idealized case 
in which V = (except in step eight) and J = in the low coupling regime with the same 
output when V ~ 1CT 3 in the high coupling regime, V ~ 1CT 2 in the low coupling regime and 
J ~ 1CT 2 in the low coupling regime. The maximum N that can be tolerated is taken to be 
the first iV for which this fidelity falls below 0.99. In this simulation we have taken t u = Is 
and 9 = tt/100. Figure [2] shows how the fidelities decrease as N increases. We see that by 
our definition the maximum number of atoms the system can tolerate is approximately 60. 
Squeezed states with larger numbers of atoms have been demonstrated experimentally [10] 
and as such interactions are likely to be one of the main limiting factors of this scheme. 
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FIG. 2: The fidelity between the output of scheme 2 in the idealized case (where V = 0, and J = 
in the low coupling regime) with the output in the non-idealized case (where V ~ 10~ 3 in the high 
coupling regime, V ~ 10~ 2 in the low coupling regime and J ~ 10 -2 in the low coupling regime) 
for different numbers of input atoms. Here 6 = 7r/100 and t u = Is. 

D. Metastability 

Our gyroscope scheme relies on the two counterpropagating superfuid states being stable 
at least over the timescale of the measurement that we want to make. There has been a 
lot of discussion and analysis of what conditions need to be met so that a superfluid cur- 
rent persists (or is metastable) in BECs. In a toroidal trap, for example, the condition 
for metastability is that the mean interaction energy per particle is greater than the single 
particle quantization energy h 2 /mL 2 [32]. This can be shown to be equivalent to the condi- 
tion that the superfluid velocity is less than the velocity of sound, which is just the Landau 
criterion for the metastability of superfluid flow. Persistent superfluid flows have recently 
been experimentally observed for BECs in toroidal traps [43j. 

There has also been an analysis of the metastability requirements for a system that is 
directly relevant to our scheme, i.e. a superfluid Bose gas trapped in a ring optical lattice 
[44] . In this work, it was shown that any disorder due to energy offsets, e, on the lattice 
sites in the ring can cause dissipation by providing a coupling pathway between states 
with opposite quasi-momentum. If, however, there are interactions between the atoms, the 
energy levels are shifted and this mismatch means that the different supercurrent states 
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become effectively decoupled and so the supercurrent should persist. Kolovsky showed that 
for states with low values of quasimomentum (such as the ±1 unit of angular momentum 
states that we consider) the minimum strength of the interactions must be V m [ Q ~ 8e/N. 
Taking parameters from our scheme, in the strong tunneling region, we have J ~ 10Hz and 
V ~ 10~ 3 Hz. This means that we require e < NV/8 ~ 10 _2 Hz, for N « 60. By stabilizing 
the energy offsets to this level (or increasing the interactions) persistent supercurrents should 
be able to be achieved. 

E. Comparison with other schemes 

At this point we note that our precision analysis is for the case of a single shot, i.e. N 
atoms are loaded into the lattice and a single measurement is made over time t w correspond- 
ing to a particle flux of N/t u . In reality, the results of many runs will be combined to give 
a measurement of the rotation. Suppose we repeat the measurement n times to give a total 
integration time of r = nt^jlTj. In this case, we get 



Substituting in approximate values for our setup in the strong coupling regime (i.e. J ~ 
10Hz, N « 60, t u « Is and L = 2tt x 20/wn [22]), we get S « 10~ 3 rads -1 /VHz. This 
compares unfavorably with other atom interferometry schemes which can achieve sensitivities 
better than 10 -8 rads _1 /vTIz [45J. These other schemes achieve improved sensitivities by 
having much larger particle fluxes (e.g. 6 x 10 8 atoms/s) and much larger areas enclosed by 
their interferometer paths (e.g. 22 mm 2 ) |45j. 

The scheme presented here is therefore unlikely to challenge the overall precision offered 
by other techniques, except perhaps in specialised cases where the number of atoms available 
is restricted to a small number or the area of the interferometer must be very small. The 
main interest of this scheme, however, is that it proposes a means of creating macroscopic 
superpositions of persistent superfluid flows and that these could show evidence of Heisenberg 
scaling of measurement precision. This in itself would be of fundamental interest. In order 




(35) 



where the short-time sensitivity is given by, 




(36) 
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to improve its short-term sensitivity, however, it is likely to be difficult to create entangled 
states with very large numbers of particles, so a different configuration would need to be 
used to greatly enhance the enclosed area of the interferometer. 

Alternatively nonlinear couplings of the atoms and the rotation could potentially be used 
to achieve precision scalings of N~ 3 ^ 2 using unentangled atoms (T9J EH]. The reliance of 
such a scheme on unentangled atoms would allow for much larger atom numbers than in 
our scheme and as such seems to offer a possible way of improving upon current gyroscope 
precisions. However, the challenge of this scheme would be coupling the angular momentum 
and the atomic scattering lengths. 

VII. CONCLUSION 

We have presented three schemes to measure small rotations applied to a ring of lattice 
sites by creating superpositions of ultra-cold atoms flowing in opposite directions around 
the ring. Two of these schemes are capable of Heisenberg limited precision measurements 
where the precision scales as 1/JV. The two schemes use different entangled states. While 
the scheme that used a NOON state gave slightly better precision in the idealized case, after 
consideration of experimental limitations it was shown that the bat state is likely to be the 
preferred candidate largely due to its robustness to particle loss. Importantly the bat state 
outperformed the case of unentangled particles for modest loss rates. The effects of non-zero 
interactions were shown to limit the preferred scheme to approximately 60 atoms and as such 
this scheme is not capable of outperforming the precision of existing atomic gyroscopes at 
present. However, the interesting result is the Heisenberg scaling of the precision. All the 
steps in this scheme should be within reach of current technologies which is promising for 
its experimental implementation. 

This work was supported by the United Kingdom EPSRC through an Advanced Re- 
search Fellowship GR/S99297/01, an RCUK Fellowship, and the EuroQUASAR programme 
EP/G028427/1. 
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